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ABSTRACT 

We show how correlated steps introduces significant contributions to the modification 
of the halo mass function in modified gravity models, taking the chameleon models as 
an example, in the framework of the excursion set approach. This correction applies to 
both Lagrangian and Eulcrian environments discussed in previous studies. Correlated 
steps also enhances the modifications due to the fifth force in the conditional mass 
function as well as the halo bias. We found that abundance and clustering measure- 
ments from different environments can provide strong constraints on the chameleon 
models. 
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1 INTRODUCTION 

The discovery of accelerated expansion of the Universe 
sparks surge of research on the possibility of modified gravity 
model (see, for example, Jain & Khoury 2010; Clifton et al. 
2012, for review). The main goal of such modification is to 
alter the large scale behaviour to explain the weakening of 
gravity - however any modifications in the gravity model 
must at the same time satisfy the tight constraints from 
the solar system test. One way to fullfill this requirement is 
to include some kinds of screening mechanisms to suppress 
the modification in local environment (high density regime) . 
In this work we focus on one particular example - one that 
modifies gravity by introducing a dynamical scale field which 
mediates a fifth force. This scalar field has the nice property 
such that the fifth force is strongly suppressed in regions 
of high matter density and hence it passes the solar system 
test. 

The chameleon model of Khoury & Weltman 
(2004); Mota & Shaw (2007) is a representative exam- 
ple (other noticeable examples are the environmentally 
dependent dilation Brax et al. (2010) and symmetron 
(Hinterbichler & Khoury 2010) models). The background 
evolution and the linear perturbations on large scale can 
be indistinguishable from the standard ACDM cosmology 
(Hu & Sawicki 2007; Li & Barrow 2007; Brax et al. 2008; 
Li & Zhao 2009); solar system test is satisfied by construc- 
tion. Hence nonlinear structure formation is the only regime 
where effects of such models would possibly be detected. 
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A number of studies (Oyaizu 2008; Oyaizu et al. 2008; 
Schmidt et al. 2009; Li & Zhao 2009, 2010; Li & Barrow 
2011; Zhao, Li & Koyama 2011; Li et al. 2011; Brax et al. 
2011; Davis et al. 2011) employed N-body numerical sim- 
ulations to study nonlinear structure formation - however 
high resolution simulations with cosmological volume are 
still challenging due to the non linear equation governing 
the scalar field. 

This work aims at investigating the effect of chameleon 
mechanism on large scale structure. We apply the excur- 
sion set approach (Bond et al. 1991; Mo & White 1996; 
Sheth & Tormen 1999) to compute the halo mass function in 
the same line as previous analyses done by Li & Efstathiou 
(2012); Li & Lam (2012). Li & Efstathiou (2012) first illus- 
trated the idea of extending the excursion set approach 
to calculate halo mass function in chameleon models (see 
Brax, Rosenfeld & Sterr (2010) for an earlier work in this 
direction, where the authors studied the spherical collapse 
in chameleon models in detail). They did so by assum- 
ing a fixed Lagrangian environment. Li & Lam (2012) im- 
proved this approach by introducing an Eulerian environ- 
ment. This choice avoids the unphysical requirement that 
all environments having the same mass (which also places 
an upper limit of the halo mass). The calculation involves 
two first crossing distributions, one for the Eulerian environ- 
ment and the other one for the (modified) halo formation 
barrier. Both studies compute the effect of the chameleon 
mechanism assuming uncorrrelated steps in the excursion 
set approach - it is the main focus of the first half of the 
current work to investigate how correlated steps in the ex- 
cursion set approach interact with the fifth force and the 
chameleon mechanism. Recently there are renewed inter- 
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ests in going beyond the uncorrected steps excursion set 
in GR+ACDM cosmology to compute the halo/void abun- 
dance (Maggiore & Riotto 2010; Paranjape, Lam & Sheth 
2012a,b; Musso & Sheth 2012) and the halo bias (Ma et al. 
2011; Paranjape & Sheth 2012; Musso & Sheth 2012). 

In addition to the unconditional mass function, the ex- 
cursion set approach lays down the framework to calculate 
the conditional mass function as well as the halo bias. The 
latter is of particular importance since it relates observables 
(distribution of halos) to the underlying matter distribution. 
Since the mass function in chameleon models depends on the 
environment density, one may expect the conditional mass 
function (and hence the halo bias) to be modified. The sec- 
ond part of this paper discusses how the conditional mass 
function as well as the halo bias depend on the chameleon 
model. 

This paper is organised as follows: in § 2 we briefly re- 
view the theoretical model to be considered and summarise 
its main ingredients. The effect of the chameleon mecha- 
nism on the halo mass function within the framework of 
excursion set approach is presented in § 3, where we also in- 
clude analytical approximation for correlated steps in § 3.4. 
§ 4 discusses the conditional mass function and halo bias in 
chameleon models where we show measurements in under- 
dense environment provide a good test for chameleon signa- 
tures. Finally we conclude in § 5. 



2 THE CHAMELEON THEORY 

This section lays down the theoretical framework for inves- 
tigating the effects of coupled scalar field(s) in cosmology. 
We shall present the relevant general field equations in § 2.1, 
and then specify the models analysed in this paper in § 2.2. 



2.1 Cosmology with a Coupled Scalar Field 

The equations presented in this sub-section can be found 
in (Li & Zhao 2009, 2010; Li & Barrow 2011), and are pre- 
sented here only to make this work self-contained. 
We start from a Lagrangian density 



R 

V a (pV a tp 



+ V{ip)-C(ip)£uM+£-s, 



(1) 



in which R is the Ricci scalar, k = 8irG with G being the 
gravitational constant, £dm and £s are respectively the 
Lagrangian densities for dark matter and standard model 
fields, ip is the scalar field and V(ip) its potential; the cou- 
pling function C(ip) characterises the coupling between ip 
and dark matter. Given the functional forms for V(p) and 
C(ip) a coupled scalar field model is then fully specified. 

Varying the total action with respect to the metric g a b, 
we obtain the following expression for the total energy mo- 
mentum tensor in this model: 



T a b 



VaipVtlfi — gab 



-X7 c X7 c p-V(ip) 



+C(<p)TT + TL 



(2) 

where T^ M and are the energy momentum tensors for 
(uncoupled) dark matter and standard model fields. The 
existence of the scalar field and its coupling change the form 



of the energy momentum tensor leading to potential changes 
in the background cosmology and structure formation. 

The coupling to a scalar field produces a direct interac- 
tion (fifth force) between dark matter particles due to the 
exchange of scalar quanta. This is best illustrated by the 
geodesic equation for dark matter particles 



df 2 



= -V<(>- 



dp) 



V<p, 



(3) 



where r is the position vector, t the (physical) time, <j> the 
Newtonian potential and V is the spatial derivative. C v = 
dC/dip. The second term in the right hand side is the fifth 
force and only exists for coupled matter species (dark matter 
in our model). The fifth force also changes the clustering 
properties of the dark matter. 

To solve the above two equations we need to know both 
the time evolution and the spatial distribution of ip, i.e. 
we need the solutions to the scalar field equation of motion 
(EOM) 

dV(ip) dC(yQ 

V V a ip-\ : h PDM : = 0, 

dip (lip 

or equivalently 

dip 

where we have defined 
V e ff(ip) = V{ip) + pdmC(ip). 



(4) 



(5) 



(6) 

The background evolution of ip can be solved easily given 
the present day value of pdm since pdm oc a -3 . We can 
then divide ip into two parts, ip = ip + Sip, where tp is the 
background value and Sip is its (not necessarily small nor 
linear) perturbation, and subtract the background part of 
the scalar field equation of motion from the full equation 
to obtain the equation of motion for Sip. In the quasi-static 
limit in which we can neglect time derivatives of Sip as com- 
pared with its spatial derivatives (which turns out to be a 
good approximation on galactic and cluster scales), we find 



^2 dCjip) 

V ip = — ; pdm 

dip 



dc(v0 dvfr) 

~Pdm + 



dip 



dp 



dV(y>) ; 
dip 



(7) 



where pdm is the background dark matter density. 

The computation of the scalar field ip using the above 
equation then completes the computation of the source term 
for the Poisson equation 



K r 

= 2^ 



+ 3p t , 



[C(ip)puM + pb- 2V(ip)} 



(8) 



where pe is the baryon density (we have neglected the ki- 
netic energy of the scalar field because it is always very small 
for the model studied here). 



2.2 Specification of Model 

As mentioned above, to fully fix a model we need to specify 
the functional forms of V(ip) and C(ip). Here we will use the 
models investigated by Li & Zhao (2009, 2010); Li (2011), 
with 



C(ip) = exp(7^<p), 



(9) 



© 2012 RAS, MNRAS 000, 1-12 



Excursion set with correlated steps for modified gravity 3 



and 

V(<p) = tz f „ ■ (10) 

[1 - exp(-V«¥>)] 

In the above A is a parameter of mass dimension four and 
is of order the present dark energy density (<p plays the role 
of dark energy in the models). 7,0 are dimensionless pa- 
rameters controlling the strength of the coupling and the 
steepness of the potentials respectively. 

We shall choose a « 1 and 7 > as in Li & Zhao 
(2009, 2010), ensuring that V e ff has a global minimum 
close to tp = and d 2 V e ff(ip)/dip 2 = mt, at this mini- 
mum is very large in high density regions. There are two 
consequences of these choices of model parameters: (1) ip 
is trapped close to zero throughout cosmic history so that 
V(ip) ~ A behaves as a cosmological constant; (2) the fifth 
force is strongly suppressed in high density regions where <p 
acquires a large mass, m 2 , 3> H 2 (H is the Hubble expan- 
sion rate), and thus the fifth force cannot propagate far. The 
suppression of the fifth force is even stronger at early times, 
and thus its influence on structure formation occurs mainly 
at late times. The environment-dependent behaviour of the 
scalar field was first investigated by Khoury & Weltman 
(2004); Mota & Shaw (2007), and is often referred to as the 
'chameleon effect'. 
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Figure 1. Upper panel: Halo formation barrier for various <5cnv 
from -1.2 to 1.2 with increment 0.2. Constant barrier at the top is 
the case with no fifth force (i.e., ACDM). Lower panel: The height 
need to climb to reach the halo formation barrier for various Ocnv ■ 



3 FIRST CROSSING PROBABILITY WITH 
CORRELATED STEPS IN CHAMELEON 
MODELS 

3.1 Terminology 

In the excursion set approach the calculation of the halo 
mass function dn /dm is mapped to the computation of the 
first crossing distribution f(S) across some prescribed bar- 
riers where 

™2ldm = f(S)dS. (11) 
p dm J w v ' 

In the following we use the two terms interchangeably. The 
variance of the matter fluctuation field smoothed on scale R 
is given by 

where W is the smoothing window function and P(k) is the 
matter power spectrum linear extrapolated to the present 
time. In hierarchical models S is a monotonic decreasing 
function of R and the smoothing scale in Lagrangian space 
R, the total mass enclosed within this scale M, and the 
variance S are equivalent quantities and can be used inter- 
changeably. 

3.2 Background 

Recent work (Li & Efstathiou 2012; Li & Lam 2012) demon- 
strated how the chameleon-type fifth force modifies the first 
crossing probability and the associated halo mass function 
for Lagrangian and Eulerian environment. This modification 
is due to the dependence of the halo formation barrier height 
S C (S) on the matter density of its surrounding environment, 



^onv 1 - Since the fifth force is strongest when the environ- 
ment density is low, the barrier for halo formation is lower 
where <5 cnv is small. The upper panel in figure 1 shows the 
halo formation barrier as a function of S for various values 
of S env assuming that the spherical collapse model describes 
halo formation, and confirms that the fifth force is stronger 
in underdense environment where the associated barrier is 
lower. The actual height that the random walks with various 
S C nv need to climb to reach the modified barrier is shown in 
the lower panel - although the barrier in underdense envi- 
ronment is lower, the barrier to overcome is actually higher. 

Li & Efstathiou (2012); Li & Lam (2012), when apply- 
ing the excursion set approach in chameleon models, focused 
on the case where the excursion set is performed with uncor- 
rected steps with two different definitions of environment. 
In both cases the modifications in the first crossing distribu- 
tion in chameleon models are solely due to the change in the 
barrier. The situation is more complicated when the steps in 
the random walk are correlated: two additional effects can 
modify the first crossing distribution. Firstly, the distribu- 
tion of 8 env for Eulerian environments is modified and, since 
the halo formation barrier depends on 5 env , the first crossing 
probability will be modified as well. Note however that when 
the environment is defined as Lagrangian, the distribution 
of <5 e nv is unchanged. Secondly, there will be non-trivial cor- 
relations between the environment density contrast 5 env and 
the density contrast S at different S. 

To investigate how these two effects modify the first 

1 Note that in this paper we use <5 e nv to denote the initial density 
perturbation in the environment linearly extrapolated to today 
assuming a ACDM model. Knowing <5 e nv it is straightforward to 
calculate the true environment density in subsequent times using 
the assumed evolution model. 
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Figure 2. (Colour Online) Distribution of environmental den- 
sity for Eulerian (top panel) and Lagrangian environments (lower 
panel). The solid, dashed and dotted curves are results of tophat, 
Gaussian and sharp- k window functions respectively. 



crossing distribution in models with chameleon mechanism, 
we use monte-carlo simulations to study the first crossing 
distribution for three choices of smoothing window func- 
tions following the procedures in Bond et al. (1991). Our 
convention is such that the normalisation of the power spec- 
trum is set by demanding as — 0.81 for the tophat filter. 
The sharp-A smoothing window function generates random 
walks with uncorrelated steps and it is common practice to 
use the S — M mapping of the tophat window function to 
relate the variance and the smoothing scale for the sharp- 
Ac smoothing window function. We consider ACDM power 
spectrum and two environment scales: an Eulerian environ- 
ment, which denoted by (, is the environment at late time; 
and a Lagrangian environment, which denoted by £, is the 
environment in the initial conditions. We assume both en- 
vironments are spherical regions that share common cen- 
ters as the proto-halos. The actual scale of the environment 
should be approximately the Compton wavelength of the 
scalar field at late time. Numerical simulation results sug- 
gest that an Eulerian scale of £ = 5 ~ 8/i -1 Mpc should be 
chosen (Li, Zhao & Koyama 2012). In this work we consider 
an Eulerian evironment £ = 5/i~ 1 Mpc and a Lagrangian en- 
vironment £ = 8/i~ 1 Mpc, and set (7, a) = (0.5, 10 -6 ) as the 
chameleon model parameters. In the language of excursion 
set formalism, the Lagrangian environment barrier (fixing 
S) corresponds to a vertical barrier in the 5i — S plane while 
the Eulerian environment barrier is given by the spherical 
collapse model (Bernardeau 1994; Sheth 1998): 



6eui(5') = S c o 



-i/<W 



(13) 



where p is the background density. 

Figure 2 shows the distribution of <5 en v for £ = 5 Mpc//i 
(upper panel) and £ = 8 Mpc/h (lower panel). Three his- 
tograms corresponding to tophat (solid) , Gaussian (dashed) , 



and sharp-A (dotted) window functions are shown in each 
panel. We ran a sample of random walk whose correla- 
tion satisfies the corresponding smoothing filter function. 
We then record the height of the random walk at which 
it first crosses the environmental barrier. In the case of 
Lagrangian environment we set S anv = S c o if the random 
walk first crosses S c o on scale bigger than £. Since we use 
the same S — M mapping for the tophat window function 
and the sharp-A; window function, choosing a Lagrangian 
scale fixes S and hence the two will have the same distribu- 
tion. The difference between the distributions of Gaussian 
smoothing and the others in the case of Lagrangian envi- 
ronment is due to our choice of power spectrum normal- 
ization - for £ = 8 Mpc//i, lg(f) = 1.1 for the Gaussian 
filter while lg(^) = 0.62 for the tophat and sharp-A; filters. 
Here v = S^ /S. Since the Lagrangian 5 Bnv distribution is 
the same for the tophat and sharp-A; filters, any difference 
in their first crossing distributions due to the fifth force has 
to come from the non-trivial correlations between the height 
of the random walk and S onv (see below). On the other hand, 
the distributions in Eulerian environment are different for all 
three smoothing windows: the sharp-A; filter peaks at high- 
est S e nv while the Gaussian filter peaks at lowest <5 onv . As 
a result the effect of the fifth force is stronger when the 
Gaussian window function is used - the halo formation bar- 
rier for walks with Gaussian filter is lower than that of un- 
correlated walks. Analytical approximations to evaluate the 
first crossing of Eulerian environment are available for both 
uncorrelated (Zhang & Hui 2006; Lam & Sheth 2009) and 
correlated steps (Musso & Sheth 2012). 

3.3 Monte-Carlo simulations 

The top panel in figure 3 shows the first crossing distribu- 
tions in chameleon models for both correlated and uncor- 
related walks. The histograms show the distribution when 
there is no fifth force: it is the first crossing for a constant 
barrier at S c = 1.676. The line labels are the same as in fig- 
ure 2. Symbols represent the first crossing distributions for 
different combinations of smoothing window functions and 
choices of environment, as indicated in the legends. The first 
crossing distributions change significantly when the random 
walks switch from uncorrelated steps (sharp-A) to the cor- 
related one (Gaussian and tophat); different definitions of 
the environment have subdominant effect. To better visu- 
alise the effect of different choices of the environment, the 
lower panel shows the relative difference in the distribution 
to the corresponding constant barrier distribution. 

The modifications in the first crossing distribution due 
to the chameleon-type fifth force are quite different for ran- 
dom walks with correlated steps (both Gaussian and tophat 
filters) versus walks with uncorrelated steps. For example 
the modifications of the correlated steps are displaced from 
those of the uncorrelated steps for \g(u) < 0.1. Since the 
modification in the halo mass function is around 20%, this 
~ 5% displacement is a big contribution. Consider for exam- 
ple the tophat (solid squares) and the sharp-A; (open squares) 
Lagrangian environments - since they have the same distri- 
bution p(5 C nv), the difference in the predicted modification 
in the first crossing distributions must be due to the nontriv- 
ial correlations between different smoothing scales when the 
window function is tophat. One can imagine the following 
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Figure 3. (Colour Online) Comparison of the first crossing distri- 
bution in chameleon models for correlated as well as uncorrclatcd 
steps. Three sets of histograms are the first crossing probability 
for the ACDM model, a constant barrier S c = 1.676, for different 
smoothing window functions: from top to bottom are sharp-fc fil- 
ter (uncorrelated walks), tophat filter, and Gaussian filter. The 
symbols are the first crossing distribution when there is the fifth 
force for £ = 5h~ 1 VLpc and £ = 8h _1 Mpc. See the legends for 
the details. The bottom panel shows the fractional change in the 
first crossing probability with respect to the corresponding con- 
stant barrier first crossing probability. Only error bars for red 
solid triangles are included for clarity. 



scenario: random walks with correlated steps are less likely 
to have dramatic fluctuations compared to those with uncor- 
related steps - when S cnv < 0, the enhanced fifth force lowers 
the barrier <5 C - in both uncorrelated and correlated cases, 
it is easier to cross the barrier. However correlated walks 
are not likely have sudden jumps, hence the modification 
in the first crossing comes only at bigger S: it explains the 
difference between the solid (tophat Lagrangian) and open 
(sharp-fc Lagrangian) symbols for lg(z^) > 0.3 (Lagrangian 
scale £ = 8 Mpc/h corresponds to lg(f ) = 0.62). For bigger S 
the correlated walks do not need have sudden jumps to cross 
the barrier - hence the increase in the first crossing proba- 
bility. On the other hand, the increase for the uncorrelated 
walks is smaller since some of those walks already cross the 
barrier at smaller S. For <5 cnv > 0, the fifth force is weak and 
only slightly enhances structure formation. In this case walks 
with correlated steps may upcross the barrier sooner than 
walks with uncorrelated steps (since uncorrelated walks can 
have sudden decrease in height). It will result in higher prob- 
ability in first crossing distribution for correlated steps for 
S near the chameleon environment - however this modifica- 
tion in the first crossing distribution relative to the GR case 
is small due to the weakening of the fifth force. Recall that it 
is equally likely to have overdense and underdense environ- 
ment in Lagrangian environment, the net effect is correlated 



steps show a smaller change in first crossing probability for 
small s but bigger change for big s. 

The change in first crossing distribution for Gaussian 
filter and Lagrangian environment is strong in high v. It is 
consistent with the change in the distribution of S cnv shown 
in figure 2: since p(S cnv ) has a narrower distribution with a 
peak at <5 cnv = 0, the overdense environments in the Gaus- 
sian case generally have lower S cnv and so the fifth force is 
stronger (lower barrier) compared to the walks with tophat 
or sharp-fc filter. On the other hand the change at low mass 
end is similar to others - although there are more walks that 
have underdense Lagrangian environments for the other two 
filters, those walks need to climb a huge distance to cross 
their respective modified barriers and do not make a signif- 
icant contribution in the first crossing distribution. 

3.4 Analytical approximation 

In this subsection we apply the analytical approximation 
proposed by Musso & Sheth (2012) to estimate the first 
crossing distribution in chameleon models. In this analytical 
approximation the first crossing probability at s is approxi- 
mated by both the height and its rate of change at s — As. 
The first crossing distribution across barrier b(s) 2 is 

dS' p(S'\b,s)(S' -b'), (14) 

where the dash denotes derivative with respect to s, p(<5'|&, s) 
is the conditional probability density of the rate of change 
in height 5' at s given that the walk's height is b at the same 
scale s, and p(8, s) is the probability density that the walk 
has height b at the smoothing scale s. 

When the Lagrangian environment is used in the 
chameleon models, the first crossing distribution of the mod- 
ified barrier b — S c (s, S env ) at s is 

dS cnv p(8 eBV ,£)p(b,s\5env,0 

-oo 

/oo 
dS' p{S'\b,s;5 cnv ,Z)(5' -b'), (15) 

where p(5 env ,£) is the distribution of environmental density 
contrast at scale £. Strictly speaking p(<5 onv ,£) should be 
replaced by the condition probability that S never crosses 
5 c o for scales bigger than £, but the scales we apply the 
Lagrangian environment is big enough so that the difference 
between the two is small. 

Figure 4 shows the comparison of the analytical approx- 
imation and the monte-carlo results for Lagrangian environ- 
ment of £ = 8Mpc//i. The upper panel shows the results of 
Gaussian window function while the lower panel shows the 
tophat window function. In each panel the top half shows the 
first crossing distribution from analytical prediction (curve) 
and the monte-carlo result (symbols), and the bottom half 
shows the modifications to the first crossing distribution 
relative to GR. The analytical approximation matches the 
monte-carlo simulation for lg(i^) > - for smaller v, multiple 
crossings are more frequent and the approximation breaks 
down. 

2 In this subsection we use upper and lower letters to denote 
different smoothing scales and random walk heights. 
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Figure 4. Analytical approximation for first crossing distribu- 
tion in chameleon models for correlated steps. The upper and 
the lower figure show the results using Lagrangian environment 
£ = 8 Mpc/h for Gaussian and tophat window filters respectively. 
In each figure the curve is the analytical approximation while the 
data points are monte-carlo simulation results. The upper half of 
each figure shows the first crossing distribution and the lower half 
shows the ratio to the constant barrier distribution. The vertical 
dotted lines show the scales of the Lagrangian environment. 

The first crossing distribution for Eulerian environment 
contains two first crossings: first to cross the environment 
barrier, then the halo formation barrier. Applying the ap- 
proximation to both crossings yields 

PS poo 

/ Eul ( S ,C) = / dS P (B,S) dA'p(A'|B,5)(A'-B') 

JO J B' 

xp(b,s\B,S; A') 

r oo 

x J dS' P (5'\b,s;B,S;A')(S'-b'), (16) 

where we denote the Eulerian environment barrier as B(S) 
and A' is the derivative of 5 em at S. The dependence of 
the halo barrier b(s) on B(S) is implicit. The first integral 
includes walks with different Eulerian environment density 
contrast (up to s) while the second integral corresponds to 
the first crossing of the Eulerian barrier at S. Here one 
should have replaced the upper limit of the integration such 
that the walk would not cross S c o at 5* - in practice the 
probability of having a big jump is not likely so using the 
above approximation does not alter the result. Finally the 
rest corresponds to the first crossing of the modified halo 
barrier at s given the walk first crossed the Eulerian barrier 
at S with an increase in height A' when going from S — AS 
to S. The results are shown in figure 5. The analytical ap- 
proximation works relatively well for big v (> 1) but breaks 
down for small u, as in the case of Lagrangian environment. 

The conditional probability distributions in this subsec- 
tion are also Gaussian distributed with mean and variance 



Figure 5. Similar to figure 4 but using Eulerian environment 
C = 5 Mpc/h. 

given by 

71 

E(x\yi,y 2 ,--- ,y n ) = ^(xyijA^yj (17) 

71 

Vax(x\yi,y 2 , ■ ■ • , y„) = Var(x) - ^(xy 3 ) 2 A^ 

3=1 

-2j2{xyj)(xy k )Aj k -, (is) 

where A is the covariance matrix of and we assume 
all variables have unconditional means vanish in the above 
formuale. 



4 CONDITIONAL MASS FUNCTION AND 
HALO BIAS IN CHAMELEON MODELS 

In this section we will describe the effect of the chameleon- 
type fifth force on the conditional mass function and the 
halo bias. As explained in the previous section, the fifth 
force is stronger when the environment is less dense - ef- 
fectively the halo formation criterion is lower and hence the 
first crossing of halo formation barrier is easier. One may 
then ask if this effect will be more significant when we look 
at the conditional mass function, especially for underdense 
regions. Here we will use Monte Carlo simulations to inves- 
tigate how the conditional mass function responses to the 
fifth force introduced by the chameleon models. 

We now have two environments: one is used to solve 
the fifth force in the chameleon model (we will call it 
the chameleon environment hereafter and we choose £ = 
5/i _1 Mpc and £ = 8/i _1 Mpc) and the other is a large-scale 
environment S = So on which the condition is defined. We 
will choose So to be small (hence very large scale) so that 
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Figure 6. (Colour Online) Comparison of the conditional first crossing probability obtained from Monte Carlo simulations for uncor- 
related and correlated steps by selecting subset of walks that have various values of <5o at So = 0.15 2 <5 2 . In each panel the histograms 
indicate the first crossing probability for the constant barrier case (ACDM) while the the solid symbols and the open symbols represent 
chameleon models with Q = 5h~ x Mpc and £ = 8h~ 1 Mpc respectively. Three sets of probability are included in each panel: |5o| (highest 
amplitude histogram at lg(v) = 0.5 and squares); — |<5rj| (lowest amplitude histogram at lg(f) = 0.5 and pentagons); and unconditional 
probability (the intermediate histogram and triangles). Different panels use different smoothing window function: (from top to bottom): 
uncorrelated steps; tophat window function; Gaussian window function. Only error bars for the histogram for positive <5o are included 
for clarity. 
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Figure 7. (Colour Online) Conditional first crossing probability similar to figure 6 the condition is set at S = 0.25 2 <5 2 . 



it is safe to assume that the range of the fifth force is much 
smaller than the scale corresponding to So. 



4.1 Conditional first crossing probability in 
chameleon models 



In this subsection we will look at the conditional first cross- 
ing probability for three different smoothing window func- 



© 2012 RAS, MNRAS 000, 1-12 



8 T. Y. Lam & B. Li 



tions. The conditional first crossing probability is 

f(S\8 ,So) = f dS cnv f(S\5 cnv ;So,So)P(S cnv \5 ,So), (19) 

where f(S\5 e nv', So, So) is the conditional first crossing prob- 
ability across the 5 C (S, S env ) barrier given that it first crosses 
the environment barrier B(s) at 5 cnv and it has S = So at So- 
P{S C nv\So, So) is the conditional first crossing probability of 
crossing the environment barrier at S cnv . Changing the large 
scale environment can have two effects: 

(i) It modifies the distribution of <5 cnv : when So < 0, it is 
more likely to have less dense chameleon environment (5 en v 
is smaller) and hence the fifth force will be stronger. 

(ii) It also affects the first crossing distribution of the halo 
formation barrier when the random walk is non-Markovian 
(this effect is additional to the first point where the halo for- 
mation barrier is changed due to the change of environment 
distribution) . 

One may expect the second effect to be weaker than the first 
since the correlation is weaker for bigger difference in 5. In 
particular, for sharp-A; filter where the random walks have 
uncorrelated steps, the walks are Markovian and the large 
scale environment (So, So) only modifies the distribution of 
5 cnv , that is, 

/ U ncor(5|5 cnv ; So, So) = f(S\S cnv ). (20) 

We are going to quantify the effect in this section us- 
ing monte-carlo simulations. We choose two large scales 
(So = 0.15 2 5 2 and 0.25 2 <5 2 ) and measure the conditional 
first crossing distributions using different window filters. 

Figures 6 and 7 show the conditional first crossing dis- 
tributions for various choices of So at So = 0.15 2 5 2 and 
0.25 5 C respectively. In each panel the histograms indicate 
the first crossing probability for the constant barrier case 
(ACDM), while the solid and the open symbols represent 
chameleon models with £ = 5 Mpc/h and £ = 8 Mpc/h 
respectively. Three sets of distribution are included in each 
panel: conditional probability with So > (highest ampli- 
tude histogram at \g(u) — 0.5 and squares), conditional 
probability with So < (lowest amplitude histogram at 
lg(i^) = 0.5 and pentagons), and unconditional probability 
(the intermediate histogram and triangles). Different panels 
use different smoothing window functions, from top to bot- 
tom: sharp-A;, tophat, and Gaussian filters. The conditional 
mass function differs from the unconditional mass function 
and the change depends on smoothing window function as 
well as the large scale density contrast <5o: in general condi- 
tional distributions using Gaussian or tophat window func- 
tions and having more extreme values of So show more sig- 
nificant changes compared to the respective unconditional 
distribution. 

The large scale environment modifies the distribution of 
the linear density contrast at the chameleon environment, 
which in turn modifies the first crossing distribution. The 
modification depends on the smoothing window functions 
due to the difference in correlation strength between various 
scales involved (the large scale environment, the chameleon 
environment, as well as the scale at which the halo barrier 
is first crossed). We took the ratio of the change in first 
crossing distribution for conditional walks to that of the un- 
conditional walks and the results are shown in figure 8 and 



9 for S = ±0.5 and ±0.9 respectively when So = 0.25 2 <5 2 . 
The upper and lower panels of the two figures show respec- 
tively Eulerian environment and Lagrangian environment 
for the chameleon models. In each panel the ratios for So < 
are shifted upwards by 0.1 and symbols represent different 
smoothing window kernels: open symbols for sharp-fc, solid 
symbols for tophat, and starred symbols for Gaussian. 

If the change in the first crossing probability induced 
by the chameleon-type fifth force did not depend on the 
large scale environment, this ratio should be unity (indi- 
cated by the dotted horizontal line). The ratio above (be- 
low) unity indicates that those walks experience stronger 
(weaker) chameleon effect compared to the unconditional 
walks. They also make relative bigger (smaller) contribu- 
tions to the change in the unconditional first crossing dis- 
tribution. The monte-carlo results are consistent with ex- 
pectations: for walks that have slightly overdense large scale 
environment (do > 0), the ratio is always less than unity - 
these walks are more likely to have overdense <5 en v (regard- 
less of the choice of Eulerian or the Lagrangian environment) 
and hence the strength of the fifth force is weaker. The ratios 
for walks with correlated steps (solid and starred triangles) 
have bigger deviations from unity compared to the ratio for 
walks with uncorrelated steps (open triangles), indicating 
the conditional distribution of S cnv has stronger correlations 
with So when steps are correlated (otherwise the correlation 
between So and 5 at some bigger s should result in more 
first crossing). When the value of So is more extreme, the 
deviations from unity is bigger indicating further weakening 
of the fifth force. 

On the other hand for walks that have underdense large 
scale environment (So < 0), the ratio for both uncorrelated 
and correlated walks stays above unity: the change in the 
first crossing distribution due to the fifth force is stronger 
than the unconditional walks. The ratio goes significantly 
above unity for correlated steps - indicating the abundance 
of intermediate mass halos in underdense environment would 
potentially provide a strong constraint on the chameleon- 
type fifth force. 

Figure 10 shows a similar plot but with the large scale 
environment set at So = 0.15 2 <5 2 . Comparison to figure 8 in- 
dicates that while the conditional distribution p(S en v\So, So) 
still differs from the unconditional distribution p(5 en v), the 
difference is less significant - very large environment (So = 
0.15 2 5 2 o) has weaker correlations with the chameleon envi- 
ronment compared to the previous case where So = 0.25 2 <5 2 . 

In this subsection we have studied the conditional mass 
function in chameleon models. The correlation between the 
large scale environment and the chameleon environment in- 
duces further modifications in the first crossing distribution 
- this effect is strongest (in terms of the modification of the 
conditional mass function) when the large scale environment 
is underdense. The corresponding change in the halo mass 
function for intermediate mass halo shows an extra ~30% 
boost compared to the unconditional mass function. It is 
straightforward to extend the analytical approximation in 
§ 3.4 to obtain the conditional first crossing distribution by 
including an additional condition (<5o, So) in all the probabil- 
ity distributions in equations (15) and (16). This is beyond 
the scope of this work. 
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Figure 8. (Colour Online) Comparison of how the change in first 
crossing probability due to the fifth force depends on large scale 
environment (So, So)- The ratio are taken from the conditional 
distributions of figure 7a and compared to ratio taken from the 
unconditional distributions. The upper panel shows results with 
Eulcrian environment £ = 5 Mpc//i and the lower panel is La- 
grangian environment £ = 8 Mpc//i. Two sets of points in each 
panel indicate So = 0.5 (lower set) and <5o = —0.5 (upper set, 
shifted by +0.1) respectively. Each set has three symbols rep- 
resenting different window functions: open (uncorrelated walks), 
solid (tophat), starred (Gaussian). Only error bars for uncorre- 
lated walks are included as an approximate error indicator. 



4.2 Halo bias in chameleon models 

Halo bias describes the relationship between the halo over- 
density Sh and the matter overdensity 5 where the halo over- 
density can be expanded as 



(21) 



On large scale the above expansion can be truncated in the 
first few orders. In particular the linear bias term bi is com- 
monly used to fit the ratio between the halo-matter power 
spectrum Phs and the matter power spectrum P$s (or its 
square as the ratio between halo power spectrum Phh and 
Pss). The excursion set approach provides a natural way 
to relate halo abundance and halo bias - by starting the 
random walks at some prescribed locations rather the ori- 
gin (Mo & White 1996; Sheth & Tormen 1999). The relative 
difference of this conditional first crossing probability to the 
unconditional one gives the left hand side of equation (21). 

Recently there are several studies (Ma et al. 2011; 
Paranjape & Sheth 2012) on evaluating the halo bias using 
the excursion set approach with correlated steps. Both of 
the analyses (using very different methods) found that the 
halo bias with correlated steps is stronger than the uncor- 
related case. Paranjape & Sheth (2012) also suggests that, 
when correlated steps are considered, the halo bias corn- 



Figure 9. (Colour Online) Conditional first crossing probabil- 
ity similar to figure 8 (note the change in y-axis range) but the 
subsets of walks pass through <5q = ±0.9. 



puted from the excursion set approach is different from the 
one obtained by taking the ratio of the halo-matter power 
spectrum and the matter power spectrum. While we are in- 
terested in computing the halo bias with correlated steps, 
we do not make the distinction here and we will present the 
halo bias evaluated within the excursion set framework. 

In the chameleon models the halo formation barrier in 
the excursion set approach depends on the environment den- 
sity. To estimate the halo bias we use the conditional first 
crossing distributions obtained in the previous subsection 
and compute the halo bias from equation (21) where 



b\8 c o = 



f(S\S ,S ] 



f(S) 



- 1 



~So~' 



(22) 



assuming high order terms can be neglected. We chose the 
conditional first crossing distribution where So = 0.15 2 5^. 
and So — ±0.35 and ±0.5 and the results are shown in fig- 
ure 11. The three panels show results using different smooth- 
ing window functions and the results for chameleon models 
are shifted by ±4. 

Halo bias using more correlated window functions has 
stronger dependence on the value of So, which can be seen 
by comparing the spread of bias among the three panels. In- 
cluding the chameleon-type fifth force modifies the halo bias 
and their relative changes are shown in figure 12 and 13 for 
the Eulerian and Lagrangian environments respectively. In 
both environment definitions the relative difference of the 
halo bias in chameleon models can be of factor of a few, re- 
gardless of the window functions used. Another noticeable 
signature is the difference of the mass range where the halo 
bias relative difference is most significant: when the tophat 
or Gaussian window filters are used, the change in halo bias 
for most overdense environment (triangles, So = 0.5) is most 
significant around lg(^) = —0.1 but it is shifted to lower 
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Figure 10. (Colour Online) Conditional first crossing probability 
similar to figure 8 but the subsets of walks pass through <5o = ±0.5 
at So = 0.15 2 <5 2 . 



Figure 12. (Colour Online) Relative difference in halo bias for 
chameleon models (Eulcrian environment) to that of GR. Tri- 
angles and solid lines are the relative difference for <5o = ±0.5 
respectively; squares and dotted lines for <5q = ±0.35. 
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Figure 11. (Colour Online) Halo bias computed from monte- 
carlo simulations using equation (22). We chose So = 0.15 2 <5 2 
and Sq = ±0.35 (solid) or ± 0.5 (open) for the conditional first 
crossing distribution where triangles (squares) represent positive 
(negative) <5q. 



mass lg(f) = —0.5 for underdense environment (solid lines, 
So = —0.5). This difference in mass range is getting smaller 
when So is less extreme (see the difference between squares 
and dotted lines for So = ±0.35). Combining the measure- 
ments of the halo bias in very overdense and underdense 
environments to the GR predictions (analytical formula or 
measurements from numerical simulations) would probably 
constrain the chameleon models. 



5 CONCLUSIONS 

In this paper we investigated the modifications in the mass 
function and the halo bias due to the fifth force in the 
framework of excursion set approach. Two different envi- 
ronment definitions (Eulerian environment at £ = 5 Mpc/h 
and Lagrangian environment £ = 8Mpc//i) are combined 
with three different smoothing window functions: Gaussian, 
tophat, and sharp-fc, to study how correlations between dif- 
ferent scales induced by the window functions interact with 
the fifth force. 

The halo formation barrier in the chameleon model de- 
pends on the environment density <5 en v - in an accompanied 
paper (Li & Lam 2012) we discussed the change in the mass 
function for Eulerian environment versus Lagrangian envi- 
ronment for uncorrelated steps (i.e. sharp-fc window filter) 
using the excursion set approach. In that case the difference 
in the two environments only modifies the distribution of 
<5env and subsequently the first crossing distribution. Addi- 
tional effect due to the non-trivial correlations between dif- 
ferent smoothing scales arises when the smoothing function 
is Gaussian or tophat. We used monte-carlo simulations to 
demonstrate the additional correlations indeed modify the 
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Figure 13. (Colour Online) Relative difference in halo bias for 
chameleon models (Lagrangian environment) to that of GR. 



value of do, we found that A \g(v) = 0.4 when So at 2-a level 
for So = 0.15 2 5 C Q. 

The effect of the correlations between different smooth- 
ing scales introduced by the use of different smoothing win- 
dow functions on halo abundance and halo bias has been the 
focus of some the recent studies. Its application with the ex- 
cursion set approach in the chameleon models results in fur- 
ther modifications in both halo abundance and the halo bias. 
In particular the change in the conditional mass function and 
the halo bias in underdense environment may provide poten- 
tially strong constraints on the chameleon models. Another 
test we are investigating is the strength of correlations be- 
tween the halo formation and the surrounding environment 
by measuring halo abundance and clustering under different 
environments in numerical simulations of chameleon models. 
This will be left for future study. 
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first crossing distribution and the contribution is significant, 
particularly in intermediate to low mass range. We then ap- 
plied the analytical formalism proposed in Musso & Sheth 
(2012) to compute the first crossing distribution for corre- 
lated steps. The analytical prediction matches the monte- 
carlo simulations well for high mass regime. 

Having shown the significance of correlated steps in the 
unconditional first crossing distribution, we examined the 
conditional mass function in chameleon models. When the 
steps in the excursion set are uncorrelated, the condition 
that random walks having passed through some So at a pre- 
scribed large scale So only alters the distribution of <5 env due 
to the Markovian nature of uncorrelated steps. Correlated 
steps (the random walk is non-Markovian) introduce an ad- 
ditional effect due to the non-trivial correlations between So 
and all other scales. We compared the change in the first 
crossing distribution due to the fifth force for conditional 
walks to that of unconditional walks and found that walks 
in underdense large scale environment experience a stronger 
modification compared to those in overdense large scale en- 
vironment - this effect is observed in all the three smoothing 
window functions (with different strengths) usually used in 
the literature. Hence combining the conditional halo mass 
function in underdense large scale environment with the un- 
conditional mass function can potentially provide a strong 
probe for modified gravity. 

Finally we investigate how the chameleon models mod- 
ify the halo bias derived from the excursion set approach. 
We found that the fifth force can modify the halo bias by a 
factor of a few at intermediate to low mass halos when un- 
derdense large scale environment is assumed. In addition we 
found that, when the correlated steps are used, the masses 
at which the halo bias is modified by the chameleon- type 
fifth force can be quite different for overdense and under- 
dense environments. While this difference depends on the 
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